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ABSTRACT 

We present ALMA Cycle I observations of the central kpc region of the luminous type-1 Seyfert galaxy NGC 
7469 with unprecedented high resolution (0.5" x 0.4" = 165 pc x 132 pc) at submillimeter wavelengths. Utiliz¬ 
ing the wide-bandwidth of ALMA, we simultaneously obtained HCN(4-3), HCO‘^(4-3), CS(7-6), and partially 
CO(3-2) line maps, as well as the 860 /im continuum. The region consists of the central ^ I" component 
and the surrounding starburst ring with a radius of ^ L5"-2.5". Several structures connect these components. 

Except for CO(3-2), these dense gas tracers are significantly concentrated towards the central ^ I", suggest¬ 
ing their suitability to probe the nuclear regions of galaxies. Their spatial distribution resembles well those 
of centimeter and mid-infrared continuum emissions, but it is anti-correlated with the optical one, indicating 
the existence of dust obscured star formation. The integrated intensity ratios of HCN(4-3)/HCO'^(4-3) and 
HCN(4-3)/CS(7-6) are higher at the AGN position than at the starburst ring, which is consistent to our previous 
findings (submm-HCN enhancement). However, the HCN(4-3)/HCO‘^(4-3) ratio at the AGN position of NGC 
7469 (LI I ±0.06) is almost half of the corresponding value of the low-luminosity type-1 Seyfert galaxy NGC 
1097 (2.0±0.2), despite the more than two orders of magnitude higher X-ray luminosity of NGC 7469. But 
the ratio is comparable to that of the close vicinity of the AGN of NGC 1068 (^ 1.5). Based on these results, 
we speculate that some other heating mechanisms than X-ray (e.g., mechanical heating due to AGN jet) can 
contribute significantly for shaping the chemical composition in NGC 1097. 

Keywords: galaxies: active — galaxies: individual (NGC 7469) — galaxies: ISM — ISM: molecules 


1. INTRODUCTION 

Investigation of the distribution, kinematics, and mass of 
dense molecular material in the centers of galaxies plays a 
key role in studying the evolution of galaxies because it is the 
reservoir of fuel for active galactic nuclei (AGNs), and also 
the site of massive star formation (starburst = SB). On the 
other hand, these heating sources significantly affect their sur¬ 


rounding gas in radiative and/or mechanical ways. For exam¬ 
ple, intense UV radiation from massive stars forms photodis¬ 
sociation regions (PDRs) around them, and X-ray dominated 
regions (XDRs), which can be larger in volume than PDRs, 
are formed in the vicinity of AGNs (e.g., Maloney et al. 1996; 
Hollenbach & Tielens 1999; Meijerink & Spaans 2005; Mei- 
jerink et al. 2007). Cosmic rays from frequent supernovae 
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(SNe) and the injection of mechanical energy due to SNe or 
AGN jet/outflow (mechanical heating) are also important for 
the chemical layout (e.g., Meijerink et al. 2011; Kazandjian 
et al. 2012, 2015; Matsushita et al. 2015). Thus, the dense 
molecular gas can provide us critical information on the na¬ 
ture of the underlying physical processes accompanying the 
non-thermal radiation in the centers of galaxies. 

Millimeter/submillimeter spectroscopic observations are 
indispensable for such a study since they can trace dense, cold 
molecular material with little dust extinction, which takes the 
bulk of gas mass in the nuclear regions of galaxies. This 
would contrast with observations of warm molecular gas such 
as H 2 1-0 S(l) in the near-IR (e.g., Hicks et al. 2009). 

Nevertheless, our knowledge on the circumnuclear dense 
molecular matter has still been limited, especially for AGNs. 
This is because there have been only a small number of 
high resolution (i.e., ^ 100 pc or higher scale) and sensitive 
observations of dense molecular gas even in nearby AGNs 
(e.g., Garcia-Burillo et al. 2010; Krips et al. 2011; Sani et al. 
2012; Matsushita et al. 2015), although several low resolu¬ 
tion surveys of dense gas have been conducted with single 
dish telescopes to address their global properties (e.g., Gao & 
Solomon 2004a,b). 

This situation is now being improved significantly by the 
advent of the Atacama Large Millimeter/submillimeter Array 
(ALMA). Indeed, the spatial distribution, kinematics, and the 
chemistry have been reported in greater detail for, e.g., NGC 
1068 (Garcia-Burillo et al. 2014; Viti et al. 2014; Takano et al. 
2014; Nakajima et al. 2015) , NGC 1097 (Izumi et al. 2013; 
Fathi et al. 2013; Martin et al. 2015), NGC 1433 (Combes 
et al. 2013), and NGC 1566 (Combes et al. 2014). Among 
these works, Izumi et al. (2013) (hereafter 113) suggested that 
HCN(4-3)/HCO+(4-3) and HCN(4-3)/CS(7-6) integrated in¬ 
tensity ratios seem to be higher in AGNs than in SB galax¬ 
ies (suhmm-HCN enhancement), although the number of the 
sample galaxies in 113 is few and the actual cause of the en¬ 
hancement is not clear. 

With the above things in mind, we here present the results 
of our high resolution ALMA band 7 observations of sub¬ 
millimeter dense gas tracers towards the central kpc region 
of NGC 7469. Because NGC 7469 hosts a luminous type-1 
AGN (L 2 -iokeV = 1.5 X 10^^ erg s“^; Liu et al. 2014), we can 
investigate the effects of AGN luminosity on the surrounding 
molecular gas with little uncertainty by directly comparing 
the results with those obtained in the type-1 low-luminosity 
(T 2 -iokeV = 6.9 X 10^^ erg s"!; Liu et al. 2014) AGN, NGC 
1097 (113). 

1.1. The target galaxy NGC 7469 

NGC 7469 is an active galaxy with a large scale (several 
kpc) stellar bar detected at near-IR (NIR) wavelength (Knapen 
et al. 2000). This galaxy is located at the distance of 70.1 
Mpc (z = 0.0164; Meixner et al. 1990) thus 1" = 330 pc, 
and hosts a luminous type-1 Seyfert nucleus as evidenced by 
broad Balmer emission lines (FWHM of H/3 = 4369 km s“^ 
Peterson et al. 2014) with time variability (e.g., Bonatto & 
Pastoriza 1990; Collier et al. 1998). Time variability is also 
confirmed by UV and X-ray observations (e.g., Kriss et al. 
2000; Nandra et al. 2000; Petrucci et al. 2004; Scott et al. 
2005), revealing this nucleus is indeed an AGN. NGC 7469 is 
also classified as a luminous infrared galaxy (LIRG) due to its 
high IR luminosity (Lg-iooo = ^ir = 10^^-^ Lq\ Sanders et al. 
2003). 


The type-1 AGN is surrounded by a luminous circumnu¬ 
clear SB ring with a radius of L5"-2.5" or 495-825 pc (e.g., 
Soifer et al. 2003; Diaz-Santos et al. 2007, see also Figure la), 
which accounts for two thirds of the bolometric luminosity of 
the galaxy (Genzel et al. 1995). Therefore, this galaxy pro¬ 
vides us with a nice opportunity to investigate how AGN and 
SB infiuence their surrounding molecular gas when observed 
at high resolution. 

The SB ring is prominent at various wavelengths including 
centimeter (Condon et al. 1991 ; Wilson et al. 1991 ; Orienti & 
Prieto 2010), far-IR (FIR, Papadopoulos & Allen 2000), mid- 
IR (MIR, Soifer et al. 2003; Gorjian et al. 2004), NIR (Genzel 
et al. 1995), optical-to-UV (Malkan et al. 1998; Diaz et al. 
2000; Diaz-Santos et al. 2007), and soft X-ray (Perez-Olea & 
Colina 1996). A tidal interaction with its neighbor IC 5283 
(^ 1.3' away from NGC 7469) is believed to have caused the 
powerful SB activity (e.g., Marquez & Moles 1994; Genzel 
et al. 1995). Dusty, young (< 100 Myr), and massive star 
clusters (individual typical mass is 10^“^ Mq) exist within the 
ring (e.g., Diaz-Santos et al. 2007). 

As for the central ^1" region, multi-wavelength studies in¬ 
cluding centimeter, K-band (2.2 /im), and 3.3 jim PAH feature 
indicate that the AGN is energetically dominant at these wave¬ 
lengths (Lonsdale et al. 2003; Genzel et al. 1995; Imanishi & 
Wada 2004), although one third of the K-band continuum is 
of stellar origin (Genzel et al. 1995). A core jet-like structure 
(e.g., Lonsdale et al. 2003; Alberdi et al. 2006) and ionized 
outfiows (Scott et al. 2005; Blustin et al. 2007) have also been 
observed in this AGN. Reverberation mapping revealed the 
mass of the central supermassive black hole (Mbh) to be ^ 1 
X 10^ Mq (Peterson et al. 2004, 2014). This Mbh is typical 
for QSOs at z = 5 on average, which will eventually grow up 
to Mbh = 10^ Mq at z = 0 (Shankar et al. 2009). In addition, 
the Eddington ratio of this AGN (~ 0.3; Petrucci et al. 2004) 
is also comparable to those of QSOs. Therefore the obtained 
spectrum of NGC 7469 could serve as a local template for 
high redshift QSOs as well. Other relevant properties of NGC 
7469 are summarized in Table 1 with the information shown 
above. 

Observations of CO molecules revealed the existence of a 
large amount of cold molecular gas (e.g., Sanders et al. 1988; 
Meixner et al. 1990; Israel 2009), as well as warm molecular 
gas (Rosenberg et al. 2015) in the central region of NGC 7469. 
The detailed spatial structure and the dynamics of the cold 
molecular gas were first investigated by high resolution (0.7") 
CO(2-l) observations (Davies et al. 2004). Their CO(2-l) 
map also showed a ring-like structure with a similar radius 
seen at other wavelengths, and a bright extended (~ 1") nu¬ 
clear region. Note that we hereafter call this kind of (100 pc 
scale) central molecular concentration a circumnuclear disk 
(CND) in general. Davies et al. (2004) also found a bar or 
a pair of spiral arms between the ring and the CND, which 
is only visible in the molecular emission, and not prominent 
at NIR wavelengths. This configuration is quite unusual (see 
Section 5). The total molecular gas mass inside the r < 2.5" 
(825 pc) region inferred from their observations is 2.7 x 10^ 
Mq . The spatial distributions of typical dense gas tracers such 
as HCN and HCO'^, on the other hand, have not been clear due 
to the lack of spatial resolution of previous studies. 

In this paper, we present the spatial distribution and line 
ratios of submillimeter dense gas tracers based on our high 
resolution (0.50" x 0.40", 1" = 330 pc) ALMA observations 
of the central kpc region of NGC 7469 involving HCN(4-3), 
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Table 1 

Properties of NGC 7469 


Parameter 

Value^ 

Reference^ 

RC3 morphology 

(R')SAB(rs)a 

(1) 

Position of the nucleus 


(1) 

q;j2ooo.o 

23^03^15.6^ 


(ij2000.0 

+08°52'26.4" 


Position angle [°] 

128 

(2) 

Inclination angle [°] 

45 

(2) 

Recession Velocity [km s“^] 

4925 (z=0.0164) 

(3) 

Adopted distance [Mpc] 

70.1 


Linear scale [pc arcsec“^ ] 

330 


Nuclear activity 

Seyfert 1 

(4) 

L2-iokeV [erg s“^] 

1.5 X 

(5) 

Lir [Lq] 

2.5 X 10“ 

(6) 

<SFR> (CND) [Mq yr-i kpc-^] 

50-100 

(7) 

Stellar age (CND) [Myr] 

110-190 

(7) 


Note. — The LSR velocity of the center of the galaxy is obtained from the CO(l-O) emission. We use this value to calculate the distance to NGC 7469 
and the linear scale, by assuming Hq = 71 km s“^ Mpc“^, Om = 0.27, and Oa = 0.73 cosmology. The <SFR> (CND) indicates the averaged star formation 
rate over the central ~ 1" region (i.e., circumnuclear disk = CND). The stellar age is also for the CND. (1) NASA/IPAC Extragalactic Database (NED, 
http : //ned. ipac . caltech. edu); (2) Davies et al. (2004); (3) Meixner et al. (1990); (4) Osterbrock & Martel (1993); (5) Liu et al. (2014); (6) Sanders 
et al. (2003); (7) Davies et al. (2007). 


HCO'^(4-3), CS(7-6), and partially CO(3-2) lines, and the un¬ 
derlying continuum emission. Our new ALMA observations 
improved this situation; we clearly separated the CND from 
the SB ring, and then obtained a spectrum of each compo¬ 
nent. The spatial resolution of this study (0.5") is compara¬ 
ble to that of the CO(2-l) observation (0.7") by Davies et al. 
(2004), so that we can directly compare the molecular gas dis¬ 
tributions and kinematics. Note that, however, the spatial res¬ 
olution of this study (^ 150 pc) would still be insufficient to 
exclusively probe the region where gas heating due to an AGN 
with NGC 7469-like luminosity dominates (see Appendix-A 
for details). 

We describe in Section 2 the details of our observations. 
An 860 jum continuum map and a nuclear spectral energy 
distribution are shown in Section 3. A full band 7 spectrum 
and some line ratios are presented in Section 4, although we 
will extensively compare the line ratios with those obtained 
in other galaxies and discuss the possible causes of the HCN- 
enhancement in our succeeding paper (Izumi et al. submit¬ 
ted). The spatial distributions of the submillimeter dense gas 
tracers are shown in Section 5. Finally, our main conclusions 
are summarized in Section 6. A kinematic study using our 
data is presented in Fathi et al. (2015). 

2. OBSERVATION AND DATA REDUCTION 

2.1. ALMA band 7 data of NGC 7469 

NGC 7469 was observed with ALMA on 2013 November 
3-4 with 28 antennas in C32-2 configuration, as a Cycle 1 
early science program (ID = 2012.1.00165.S, PI = T. Izumi). 
Baseline lengths range from 15.0 to 1284.3 m (corresponds 
to the uv distance from 17.4 to 1493.4 kA at 860 jim). The 
observations were conducted in a single pointing with a 18" 
field of view, which fully covered the nucleus and the SB ring 
(total ^ 4" in diameter). We set the phase tracking center 
to (aj 2000 . 0 , <5j2000.0) = (23‘’03"'15.64*, +08°52'25.80"). The 
receiver was tuned to cover the redshifted lines of HCN(4- 
3) (rest frequency z/rest = 354.505 GHz), HCO'^(4-3) (z/rest = 
356.734 GHz) in the upper side band, and CS(7-6) (z/rest = 
342.883 GHz), CO(3-2) (z/^est = 345.796 GHz) in the lower 
side band, both in the 2SB dual-polarization mode. Each 
spectral window has a bandwidth of 1.875 GHz, and two spec¬ 
tral windows were set to each sideband to achieve a total fre¬ 


quency coverage of ^ 7.5 GHz. The velocity spacing was 
originally 0.43 km s“^ (488 kHz) per channel, but 47 chan¬ 
nels were binned to improve the signal-to-noise (S/N) ratio, 
which resulted in a final velocity resolution of 20 km s ^ In 
the following, we express velocities in the optical convention. 
Note that we could only cover the CO(3-2) line with Vlsr > 
4900 km s“^ due to the limited band width for ALMA Cycle 1 
observations, the recession velocity (4925 km s“^, defined by 
CO(l-O) observations; Meixner et al. 1990), and the molecu¬ 
lar line width of NGC 7469 (typically full width at half maxi¬ 
mum = FWHM ^ 200 km s“^ at the nucleus). The bandpass, 
phase, and flux were calibrated with 3C454.3, 12257-^0743 
(this object is 2° away from NGC 7469), and 12232-^117, re¬ 
spectively. Weather conditions were good throughout the ob¬ 
servations with system temperatures of 100-200 K (usually 
100-150 K). The total on-source integration time was ^55 
min. 

The reduction and calibration of the data were done with 
CASA version 4.1 (McMullin et al. 2007) in standard man¬ 
ners, and the data was delivered to us from the East Asia 
ALMA Regional Center. All images of the line and con¬ 
tinuum emissions were reconstructed with the CASA task 
CLEAN (gain = 0.1, threshold =1.0 mJy, weighting = nat¬ 
ural). These images were further analyzed with MIRIAD 
(Sault et al. 1995). The achieved synthesized beams were 
typically 0.50" x 0.40", corresponding to 165 pc x 132 pc 
(^ 150 pc for the geometric mean) at the assumed distance of 
NGC 7469. The rms noise in the line channel maps is 0.55- 
0.58 mJy beam“^ at a velocity resolution of 20 km s“^ after a 
primary beam correction, and measured at areas free of line 
emission. However, in channels including CO(3-2) emission, 
the rms increases to ^ 2.3 mJy beam“\ probably due to re¬ 
maining sidelobes. For the continuum emission, the channels 
free of line emission were averaged, which resulted in an rms 
of 0.09 mJy beam“^ centered at z/rest = 349.7 GHz (860 yum), 
after combining both the LSB and USB data. This continuum 
emission was subtracted in the zzv-plane before making line 
maps. The conversion factors between Jy beam“^ and K are 
51.50, 51.60, 50.47, and 51.83 K (Jy beam-^)-! for CS(7-6), 
CO(3-2), HCN(4-3), and HCO+(4-3), respectively. 

We used the MIRIAD task MOMENT to make the integrated 
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intensity maps presented in Section 5, by computing the ze¬ 
roth moment of each data cube without any clipping. The 
velocity range included for the integration was ^ 4650-5250 
km s-' for CS(7-6), HCN(4-3), and HCO+(4-3), but the lower 
limit was 4900 km s“' for CO(3-2) due to our spectral set- 
ting restrictions. These velocity ranges can cover the whole 
spectral range obtained within our 18" field of view, except 
for CO(3-2). Throughout this paper, the pixel scale of ALMA 
images is set to 0.1", and the displayed errors indicate only 
statistical ones unless mentioned otherwise. We use line in¬ 
tensities corrected for the primary beam attenuation for quan¬ 
titative discussions. But this correction is not critical since 
most of the emission is within ^ 4". The adopted systematic 
error of the absolute flux calibration was 10 

2.2. Ancillary data 

We retrieved the Hubble Space Telescope (HST) ACS/HRC 
image (F550M) of NGC 7469 from the Hubble Legacy 
Archive (HLA)^. This image was calibrated by the HST 
pipeline products in a standard manner. We used the im¬ 
age only to check the positional consistency or shift between 
peaks in optical (V-band) and our ALMA band 7 images as 
shown in Section 3. 

3. 860 ijM continuum EMISSION 

Figure la shows the 860 pm continuum emission overlaid 
on the HST V-band (F550M) emission in the central 2 kpc 
region of NGC 7469. In this figure, we can see a bright cen¬ 
tral source, which corresponds to the CND as defined in Sec¬ 
tion 1.1, and some knots in the SB ring with a radius of ^ 
1.5" 500 pc). The peak position of the 860 pm contin¬ 

uum emission is ((Tj2ooo.o, fcooq.o) = (23^ 03”^ 15^62, -h08° 
52' 26.04"), which precisely coincides with that of the VLA 
8.4 GHz (3.5 cm) continuum emission (Condon et al. 1991; 
Orienti & Prieto 2010; their angular resolutions were 0.22"- 
0.25"). This precise match supports the accuracy of the as¬ 
trometry of ALMA. We hereafter treat this peak position as 
the position of the AGN of NGC 7469. The 860 pm fiux 
density there (marked as A in Figure la) is 5.19±0.09 mJy 
beam“^ The 860 pm knots in the SB ring (marked as B, C, 
and D in Figure la) coincide with those visible at centime¬ 
ter wavelengths (e.g., 8.4 GHz; Orienti & Prieto 2010), and 
the knots of dense molecular gas tracers such as CO(3-2), 
HCN(4-3), and HCO+(4-3) (Section 5). On the other hand, 
there is inconsistency between the peak positions of the 860 
pm and those of the V-band emissions. The magnitude of this 
positional shift is ^ 0.5", which is the same size as our syn¬ 
thesized beam and thus significant. Similar discrepancy was 
also reported between the PdBI CO(2-l) emission and HST J- 
band maps (Davies et al. 2004), and between centimeter and J, 
H, K-bands maps (Genzel et al. 1995), for example. Although 
one possible explanation for the discrepancy is that the abso¬ 
lute astrometry of HST is slightly off, it is rather likely that 
the stellar light is obscured by large amounts of dust traced by 
the 860 pm emission as also suggested by Diaz-Santos et al. 
(2007). 

We also investigated the origin (thermal or non-thermal) of 
the 860 pm continuum at the position-A by constructing the 
spectral energy distribution (SED) as shown in Figure lb. We 
used high angular resolution (comparable to, or higher than 

^ https://almascience.esc.org/documents-and-tools/ 
cycle-1/alma-proposers-guide 

^ http://hla.stsci.edu 


our ALMA band 7 resolution = 0.5") data of 4.9 GHz from 
Wilson et al. (1991), 8.4 GHz and 14.9 GHz from Orienti & 
Prieto (2010), 11.8 pm and 18.7 pm from Reunanen et al. 
(2010), and J, H, and K-bands from Genzel et al. (1995). The 
systematic fiux uncertainty is included. Although the data is 
sparse at the submilimeter to FIR range, the SED indicates 
that the centimeter continuum emission can be well repro¬ 
duced as synchrotron (a = -0.62±0.02; fiux oc z/^), whereas 
the fiux density of the 860 pm continuum is well beyond 
the expected value of the synchrotron emission, i.e., the 860 
pm continuum is dominated by thermal emission. Note that 
the percentage of the 860 pm continuum that is due to syn¬ 
chrotron is ^ 30% as estimated from the above scaling. In ad¬ 
dition, we used the CLEAN algorithm to deconvolve the con¬ 
tinuum emission in the LSB and USB of ALMA band 7 sepa¬ 
rately. Then, the fiux at position-A is 4.66±0.13 mJy beam“^ 
and 5.00±0.13 mJy beam“^ at 343.7 GHz and 356.0 GHz, 
respectively (errors indicate statistical ones only since they 
were observed simultaneously), suggesting a positive spectral 
index (/3 = 0±1; fiux oc but the number of data points 

is just two for this fitting). Therefore, we conclude that the 
860 pm continuum at the position-A mainly originates from 
thermal dust. 

4. BAND 7 SPECTRA, LINE RATIOS, AND DENSE MOLECULAR 
GAS MASS 

Here we show the ALMA band 7 spectra in Figure 2. These 
spectra were extracted from a single synthesized beam (150 
pc-resolution) placed at the positions A-D in Figure la, for 
comparison. The strongest line at the position-A is CO(3- 

2) , followed by HCN(4-3), HCO+(4-3), and CS(7-6). Other 
lines were not detected (i.e., < 3 a) or are blended with other 
lines, and their spectral positions are marked by black arrows 
in Figure 2. On the other hand, at the SB ring (positions B- 
D), HCO+(4-3) is brighter than HCN(4-3). The FWHM of 
the HCN(4-3) and HCO'^(4-3) emissions are 200 and 190 km 
s“\ 66 and 60 km s“\ 41 and 42 km s“\ and 50 and 50 km 
s“^ at the positions A-D, respectively (see also Figures 4 and 

5) . Thus, these values are consistent with each other, which 
suggests that they are tracing the same volume. We list the 
emission line parameters at these positions in Table 2 with 
other relevant quantities. 

4.1. ^HCN/HCO+ ^^^^HCN/CS 

By inspecting the spectra, one may notice that the HCN(4- 

3) /HCO’^(4-3) integrated intensity ratio (hereafter Rhcn/hco+) 
is higher than unity only at the position-A. To further clarify 
this, we list the values of Rhcn/hco+. and the HCN(4-3)/CS(7- 

6) integrated intensity ratio (hereafter Rhcn/cs ^ 113 also pro¬ 
posed this ratio to be an empirical discriminator of AGN and 
SB) at the positions A-D in Table 3. The integrated intensity 
of each line emission at each position is given in Table 2. They 
are calculated by the Gaussian fit to the data. The line ratios 
in Table 3 are consistent to the submm-HCN enhancement 
scenario first proposed by 113, i.e., AGNs tend to show higher 
Rhcn/hco+ and/or Rhcn/cs compared to SB galaxies, although 
the Rhcn/hco+ at the position-A is only slightly larger than 
unity and totally comparable to the Rhcn/hco+ of a typical 
early phase (e.g., Martin et al. 2009) SB galaxy NGC 253 
(1.1 ±0.3; 113). But the ratio is significantly higher than that 
of a typical late phase (e.g., Martin et al. 2009) SB galaxy, 
M82 (Rhcn/hco+ 0.4; 113). Note that a mechanical heating 
due to, e.g., frequent supernovae, would be responsible in the 
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Figure 1. (a) Continuum map towards the central 2 kpc region of NGC 7469 at 860 iim (contours) overlaid on the HST-V band (F550M) map (color; the intensity 
scale only indicates the counts). Contours are 3, 5, 10, 20, 30, 40, and 50 a, where 1 cr = 0.09 mJy beam“^. The maximum is 5.19 mJy beam“^. The central cross 
indicates the peak position of the VLA 8.4 GHz continuum (Condon et al. 1991), which is identical to that of the 860 /xm continuum. We regard this position 
as that of the AGN. The white filled ellipse indicates the synthesized beam at 860 /xm (0.50" x 0.40", PA = 68.6°). The AGN position and three bright knots 
in the SB ring are marked as A, B, C, and D, respectively, (b) The nuclear SED of the continuum emission in NGC 7469. We used 4.9 GHz data from Wilson 
et al. (1991), 8.4 GHz and 14.9 GHz data from Orienti & Prieto (2010), our ALMA band 7 data (red square), 18.7 /xm and 11.8 /xm data from Reunanen et al. 
(2010), and J, H, K-bands data from Genzel et al. (1995). The solid and dashed lines indicate and scalings with a = -0.62 and /3 = 1.5, respectively. The 
systematic flux errors are included in this plot. 
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Figure 2. ALMA band 7 spectra extracted with a single synthesized beam (150 pc-resolution) placed at positions-A, B, C, and D in Figure la. The flux of 
the spectral window containing CO(3-2) emission line is reduced by 20 times for a presentation. At the position-A (= AGN), we detected CO(3-2), HCN(4-3), 
HCO^(4-3), and CS(7-6) with > 3 a significance. Other lines were undetected (i.e., < 3 (t) or blended with other lines. Their expected locations are marked 
by black arrows in the top panel. The redshift of the nucleus assumed in this study is cz = 4925 km s“^ based on CO(l-O) observations (Meixner et al. 1990). 
The Gaussian fit to the observed CS(7-6), CO(3-2), HCN(4-3), and HCO^(4-3) emissions is represented by the red line. The fitting parameters can be found in 
Table 2. As for CO(3-2), we could not cover the whole spectral range in our observations. Thus, we use the peak flux density in the channel map (Figure 3) and 
the FWHM of HCN(4-3) emission line, to anyhow fit the CO(3-2) emission. However, we did not fit the CO(3-2) emission line at the positions A and B, since 
we could not find a clear turnover feature or blue part of the spectrum there. The continuum emission has been subtracted. Note that there is a dip at the lower 
frequency side of the CO(3-2) emission line at the position-B, which is due to a negative sidelobe (see also Figure 3). 
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central region of NGC 253 (Rosenberg et al. 2014a) for its 
relatively high Rhcn/hco +5 compared to M82. 

By investigating the 8.4 GHz radio flux obtained at 0.2" 
(Colina et al. 2001) and at 0.03", Davies et al. (2007) esti¬ 
mated that a nuclear supernova rate is ^ 0.1 yr“^ This rate 
corresponds to a star formation rate of ^ 15 Mq yr“^ assum¬ 
ing the Salpeter initial mass function (Salpeter 1955). Our 
relatively large spatial resolution (~ 150 pc) would pick up 
line fluxes emanating from the extended SB activity inside the 
CND (see also Appendix-A for the estimated size of the XDR 
of NGC 7469), which would impact the line ratios discussed 
here. We can not And any morphological and kinematic sig¬ 
natures of a jet-ISM interaction in the high resolution maps 
of H 2 and Bry emissions (Hicks et al. 2009; Muller-Sanchez 
et al. 2011), although Lonsdale et al. (2003) reported the ex¬ 
istence of a core radio jet-like structure. Therefore, we specu¬ 
late that the energetics of the CND of NGC 7469 is influenced 
by both the AGN (XDR) and the SB activity (PDR) when ob¬ 
served at the rather coarse resolution (150 pc in this work). 

In our succeeding paper (Izumi et al. submitted), we con¬ 
ducted non-local thermodynamic equilibrium (non-LTE) line 
modelings by the RADEX code (van der Tak et al. 2007). 
Assuming a one zone cloud, we widely changed gas kinetic 
temperature, molecular column density, molecular abundance 
ratio, and background temperature (we expect higher back¬ 
ground temperature in AGNs than in SB galaxies), but flxed 
H 2 gas density among the models (^h 2 = 10^ cm“^). This 
treatment was based on the assumption that there would be 
no systematic variation in ^h 2 between AGNs and SB galax¬ 
ies. Then, we calculated the resultant line intensity ratio, i.e., 
Rhcn/hco+ 01* ^HCN/cs^ ^ach model case. As a result, we 
found that both X(HCN)/X(HCO+) > 10 andX(HCN)/X(CS) 
>10 are required to reproduce RHrN/Hro+ > 1 and Rhcn/cs 
> 10 as are observed in NGC 7469 and NGC 1097. Eur- 
thermore, the line opacity of HCN(4-3) is moderately thick, 
whereas that of CS(7-6) is thin when Rrcn/cs ^10- Here, 
X(mol) indicates a molecular fractional abundance relative to 
H 2 . The above required abundance ratios are signiflcantly 
higher (at least by a factor of a few) than those necessary to re¬ 
produce the line ratios in SB galaxies. Therefore, these results 
suggest that a chemical abundance variation between AGNs 
and SB galaxies would be the key for the HCN-enhancement. 

In the following, we discuss possible mechanisms to real¬ 
ize the abundance variation in AGNs, although it is still un¬ 
clear, since (1) there are various ways to change the chemical 
compositions signiflcantly (e.g.. X-ray irradiation, mechani¬ 
cal heating, high temperature chemistry), and (2) the number 
of high resolution observational data of dense gas tracers is 
currently insufficient. 

In conventional XDR models (Lepp & Dalgarno 1996; 
Meijerink & Spaans 2005), we And reproducing a high 
X(HCNyX(HCO'^) such as > 10 is hard (usually it is less 
than unity), since a powerful ionizing capability in XDRs 
will favor to enhance ionic species such as HCO”^. Nonethe¬ 
less, according to the steady state, gas phase 1-dimensional 
XDR models intensively presented by Meijerink & Spaans 
(2005) and Meijerink et al. (2007), we can expect a high 
X(HCN)/X(HCO'^) and the resultant Rhcn/hco+ > 1 in a low- 
to-moderate gas density (^h 2 ^10^ cm“^) cloud with a low 
hydrogen column density (Ah < 10^^-^ cm“^) and a high X-ray 
flux (e.g., Ex > 10 erg s“^ cm“^). Eor larger Ah where the Ex 
is attenuated, they predicted Ehcn/hco+ < 1- Eollowing this 
scenario, we estimate the largest distance from the AGN that 


can keep Ex > 10 erg s“^ cm“^ is r ^ 18 pc and ~ 35 pc for 
NGC 1097 and NGC 7469, respectively, considering their 2- 
10 keV X-ray luminosities (L 2 -iokeV = 6.9 x 10"^^ erg s“^ and 

1.5 X 10^3 erg s”! for NGC 1097 and NGC 7469, respectively; 
Liu et al. 2014). We here assume that the X-ray emission is 
isotropic around the AGN itself. Attenuation due to obscuring 
material is considered for this estimation (Morrison & Mc- 
Cammon 1983; Appendix-B). But one critically weak point 
in this scenario is that the model predicted line intensities are 
extremely low (Meijerink et al. 2006, 2007), which is incon¬ 
sistent to the prominent HCN(4-3) and HCO‘^(4-3) emission 
lines in these galaxies. Eurthermore, Martin et al. (2015) ar¬ 
gued that HCN(1-0)/HCO‘^(1-0) integrated intensity ratio is 
decreased in NGC 1097 as we approach the very central re¬ 
gion (^ 2 at the edge of the CND^, but ^ 1.5 at the center), 
which seems to be incompatible with the above speculation of 
the largest distance for a high line ratio in NGC 1097. Note 
that, however, Martin et al. (2015) did not fully resolve the 
CND of NGC 1097, thus, their argument is a tentative one at 
this moment. 

Another possible scenario is the high temperature gas phase 
chemistry (Harada et al. 2010), which signiflcantly enhances 
X(HCN) especially at E > 300 K due to the activated forma¬ 
tion path from CN (see also Eigure 17 of 113). This mecha¬ 
nisms would be important to make line intensity of HCN(4-3) 
sufficiently detectable. Eurthermore, X(HCO'^) will be some¬ 
what decreased in these models since they will react with H 2 O 
(this is also a typical molecule in a high temperature environ¬ 
ment) to form HsO”^. This type of chemistry can be comple¬ 
mented to XDRs (Harada et al. 2013) where we can expect 
much higher gas temperature in such regions due to efficient 
X-ray heating than in SB environments. Indeed, some frac¬ 
tion of warm molecular gas exists in NGC 7469, as evidenced 
by the detections of near-IR H 2 (Hicks et al. 2009) and high-/ 
CO emission lines (Rosenberg et al. 2015). 

With these items, we here focus on the notable fact that 
the Ehcn/hco+ obtained at the position-A of NGC 7469 
(1.11 ±0.13) is only as half as that at the CND of NGC 1097 
(2.0±0.2; 113), despite the more than two orders of magnitude 
higher X-ray luminosity of NGC 7469. Note that NGC 1097 
also hosts a nuclear SB activity in the CND (^ a few Mq yr“^; 
Davies et al. 2007), which would contaminate the line fluxes 
obtained with the ^ 100 pc resolution beam (113). Eurther¬ 
more, Garcia-Burillo et al. (2014) found Ehcn/hco+ = L5 at 
the vicinity of the AGN in NGC 1068 (L 2 -iokeV = 6.9 x 10"^^ 
erg s“^; Panessa et al. 2006), whereas it increases to ~ 3 at the 
EastAVest-knots of its CND, which are ^ 100 pc away from 
the AGN (their spatial resolution was ^35 pc). These results 
indicate that higher X-ray luminosity does not simply lead to 
higher Ehcn/hco+- Erom the perspective of high temperature 
chemistry, the results again seem to be controversial, since the 
large difference in X-ray luminosity would compensate the ^ 

2.5 times larger spatial scale (in area) sampled in NGC 7469 
than in NGC 1097, leading to a higher X-ray energy deposi¬ 
tion rate per particle (this is an important parameter to deter¬ 
mine the thermal structure of XDRs e.g., Maloney et al. 1996) 
in NGC 7469, if the molecular gas distributes uniformly over 
the CND. In this case, we can expect a higher Ehcn/hco+ in 
NGC 7469 than in NGC 1097 by a naive view, but the actual 
result is totally opposite. 

Therefore, we presume that a non-radiative heating, namely 

^ An upper limit of the image convolved size is ~ 90 pc x 70 pc (113) 
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Table 2 

Parameters of detected lines at the peak positions of the 860 /xm continuum emission 


Emission 

(1) 

^rest 

[GHz] 

(2) 

^crit 

[cm-3] 

(3) 

Eu/k^ 

[K] 

(4) 

Peak flux 
[mJy beam“^ ] 

(5) 

vlsr 
[ km s-^l 
(6) 

Av 

[km s“^ ] 

(7) 

[Jy beam“^ km s“^ ] 
(8) 

1 

[Kkms“^] 

(9) 

CS(/=7-6) 

342.883 

i.gxio” 

65.8 

2.20±0.58 

4919.7±29.3 

196.5±27.2 

0.46±0.14 

23.7±7.1 





<1.74* 

- 

- 

<0.12+ 

<6.2+ 





1.87±0.58 

5059.9±60.6 

59.0±18.0 

0.12±0.05t 

6.0±2.6t 





1.98±0.58 

4982.2±47.9 

74.1±17.0 

0.16±0.06t 

8.0±3.0t 

CO(/=3-2) 

345.796 

8.4x103 

33.2 

- 

- 

- 

- 

- 





202±2** 

- 

- 

- 

- 





232±2** 

- 

- 

9.7±1.2t 

498.9±61.0t 





235±2** 

- 

- 

12.0±0.4t 

617.5±19.2t 

HCN(/=4-3) 

354.505 

8.5x10'’ 

42.5 

21.00±0.55 

4915.8±7.8 

199.6±6.7 

4.46±0.19 

225.1±9.5 





6.41±0.55 

4935.0±17.5 

65.5±5.9 

0.45±0.06 

22.6±2.8 





6.84±0.55 

5044.2±26.5 

41.0±5.0 

0.30±0.04 

15.0±2.0 





8.98±0.55 

4981.3±16.3 

50.1±1.5 

0.48±0.03 

24.2±1.7 

HCO+(/=4-3) 

356.734 

1.8x10'’ 

42.8 

19.38±0.58 

4914.1±5.7 

189.7±4.8 

3.91±0.15 

202.8±7.9 





9.09±0.58 

4936.6±13.0 

60.3±3.9 

0.58±0.05 

30.2±2.7 





14.50±0.58 

5045.2±8.9 

41.5±1.9 

0.64±0.04 

33.2±2.0 





18.24±0.58 

4984.4±6.1 

50.1±1.5 

0.97±0.04 

50.4±2.2 

860 /xm 

349.7 

- 

- 

5.19±0.09 

- 

- 

- 

- 





1.04±0.09 

- 

- 

- 

- 





1.23±0.09 

- 

- 

- 

- 





1.70±0.09 

- 

- 

- 

- 


Note. — Column 1: Full line name. Column 2: Line rest frequency. Column 3: Critical density of the line, calculated for = 100 K in the optically thin 
limit (Greve et al. 2009). Column 4: Upper level energy. Column 5: Peak flux density of the line and continuum emission determined from the Gaussian fits to 
the data (Figure 2). They are extracted with a single synthesized beam (150 pc resolution) placed at the positions-A (AGN-position), B, C, and D (listed vertically 
in this order) marked in Figure la. For the lines, errors were estimated from the adjacent emission-free channels except for CO(3-2). For the continuum, an error 
was estimated from the emission-free areas in the continuum image. *The 3 a upper limit is listed. **The peak flux density (i.e., flux density at the turnover 
feature) in the channel map is listed. Thus, these are not the results of the Gaussian fit to the data. The primary beam attenuation is corrected. Column 6, 7: LSR 
velocity and full width at half maximum (FWHM) of the line, which is calculated by the Gaussian fit. Column 8, 9: Velocity-integrated intensity in the unit of 
[Jy beam“^ km s“^] and [K km s“^] at the same positions as listed in Column 5, which is estimated by combining the Columns 5 and 7. "^The velocity centroid 
and the FWHM of HCN(4-3) at the corresponding position are used for the Gaussian fits. ^ Although the resultant integrated intensity is < 3 a, we list the value. 


Table 3 

HCN(4-3)/HCO+(4-3), HCN(4-3)/CS(7-6), and HCN(4-3)/CO(3-2) integrated intensity ratios in NGC 7469 


Position 

A 

B 

C 

D 

^HCN/HCO+ 

1.11 ±0.06 

0.75±0.11 

0.45±0.07 

0.48±0.04 

^HCN/CS 

9.50±2.87 

>3.65 

2.50±1.13 

3.03±1.15 

^HCN/CO 

- 

_* 

0.03±0.01 

0.04±0.01 


Note. — Row 1: The positions marked in Figure la, where line ratios are extracted with a single synthesized beam (150 pc-resolution). Row 2: HCN(4- 
3)/HCO^(4-3) integrated intensity ratio. Row 3: HCN(4-3)/CS(7-6) integrated intensity ratio. Row 4: HCN(4-3)/CO(3-2) integrated intensity ratio (tentative). 
*The peak intensity ratio at the position-B is ~ 0.03. The integrated intensities of these lines are shown in Table 2, although the values of CO(3-2) are tentative. 
The velocity ranges integrated over are shown in Section 2. The ratios are measured in the brightness temperature scale. 


mechanical heating caused by, e.g., AGN jet or outflow, is im¬ 
portant to realize the high /?hcn/hco+ of 2-3 (Kazandjian et al. 
2012, 2015) as is observed in NGC 1097 and in the EastAVest 
knots of the CND of NGC 1068. Indeed, a large radio jet 
exists in NGC 1068, which is supposed to cause shocks in 
the knots (e.g., Krips et al. 2011). Furthermore, a radio core 
jet has been discovered in NGC 1097 (e.g., Mezcua & Pri¬ 
eto 2014), although its intrinsic (deprojected) size is uncer¬ 
tain. However, we can not discard the possibility such as 
the line emitting region of NGC 1097 is extremely centrally- 
concentrated, which would result in higher energy deposi¬ 
tion rate per particle (i.e., more violent X-ray heating) than 
in NGC 7469 and the close vicinity of NGC 1068. Moreover, 
we again emphasize that the line fluxes of NGC 7469 would 
be contaminated by the emission from the extended SB re¬ 
gion inside the CND. Higher resolution observations of these 
molecular emissions are necessary to accurately map the dis¬ 
tribution of the line ratios across the CND of NGC 1097 and 
NGC 7469, and to search for a signature of shock-ISM inter¬ 
action. 


4.2. Other line emissions and their ratios 

Table 3 also includes HCN(4-3)/CO(3-2) integrated inten¬ 
sity ratios (^hcn/co) extracted at the positions C and D. We 
exclude the ratio at the positions A and B since we could 
not measure the accurate fluxes there. Admitting the uncer¬ 
tainty in CO(3-2) integrated intensity (we used the velocity 
centroid and the FWHM of HCN(4-3) emission line to flt the 
Gaussian profile to CO(3-2) line; see also Figure 2 and Ta¬ 
ble 2), the tentative Rucn/co ^1 the SB ring are well below 
the value obtained in Ori-KL (^ 0.32; Schilke et al. 1997). 
Since the line critical density is orders of magnitude different 
between HCN(4-3) and CO(3-2) (see Table 2), we can expect 
the line emitting region of HCN(4-3) is much smaller than 
that of CO(3-2). Hence, this ratio would be highly prone to 
be affected by the spatial resolution, which could lead to sig¬ 
nificantly lower Rucn/co NGC 7469 than in Ori-KL. This 
is well manifested in the case of NGC 253 (SB galaxy). When 
observed at ^ 200-300 pc resolution (comparable to our 150 
pc resolution for NGC 7469), Rucn/co ^ 0-03 (Knudsen et al. 
2007; Israel et al. 1995), which is quite similar to our result. 
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However, the ratio increase to ^ 0.08, i.e., increase by a fac¬ 
tor of a few, when observed at higher resolution (25-30 pc; 
Sakamoto et al. 2011). Then, matching the spatial resolution 
is important to compare line ratios of emission lines with to¬ 
tally different critical densities among objects. 

We also found that the H12 cN( 4-3)/H13CN(4-3) isotopic ra¬ 
tio at the position-A was > 11.9 in the brightness temperature 
scale. Note this is a peak intensity ratio and we used a 3 cr up¬ 
per limit for H^^CN(4-3). Then, adopting a ratio of 

> 40 observed in NGC 253 (Henkel et al. 2014; Martin et al. 
2010, ^ 340 pc resolution), the upper limit to the opacity of 
the H^^CN(4-3) is 3.5 at the position. For this calculation, 
we assume (1) the gas is in the LTE, (2) H^^CN is optically 
thick, (3) Hi^CN is optically thin, and (4) H^^CN and H^^CN 
emissions are coming from the same volume of gas. Alterna¬ 
tively, if we adopt a ratio of 100 obtained in Mrk 231 

(Henkel et al. 2014, ^ 19 kpc resolution), the upper limit of 
the opacity increases to 8 . 8 . In either case, the derived opti¬ 
cal depth is rather smaller than the values of the Galactic star 
forming regions (e.g., ^ 70 in Ori-KL; Schilke et al. 1997), 
likely reflecting the much more turbulent environments in the 
nuclear regions of galaxies. 

Note that vibrationally excited HCN emission, namely 
HCN(v 2 = 1 '^, /=4-3) at = 356.256 GHz, was not de- 
tected at the position-A, despite the fact that NGC 7469 is 
classifled as a LIRG (Sanders et al. 2003) hosting a lumi¬ 
nous AGN. Considering the above mentioned modest optical 
depth of H12 cN( 4-3), we can expect HCN(v 2 =l^^, 7=4-3) is 
completely optically thin, which makes itself hard to be de¬ 
tectable. Therefore, the non-detection of HCN(v 2 =l 7=4-3) 
emission does not necessarily mean that IR-pumping has little 
contribution to the observed line intensities, although we here 
simply treat that the HCN(4-3) emission is the consequence of 
the purely rotational transition. We need quantitative model¬ 
ing of the line excitation to further address this problem. We 
also note that the contribution of HC3N(39-38) emission to 
the HCN(4-3) emission is negligible at the position-A since 
HC 3 N emission is typically much fainter than that of HCN, 
and the line profile of the bluer HCN(4-3) is very similar to 
that of HCO'^(4-3), which does not suffer from line blending. 
We thus conclude that the contribution of HC3N(38-37) to the 
CO(3-2) emission is also negligible. 

4.3. Dense molecular gas mass 

In addition to Rhcn/hco+ ^hcn/cs^ we list line lu¬ 
minosities of HCN(4-3) and HCO‘^(4-3) extracted at the 
position-A and the whole CND (i.e., central 1" diameter re¬ 
gion) in Table 4 as well since these lines are relatively strong 
enough to be detectable even in high redshift galaxies (e.g., 
Barvainis et al. 1997 ; Riechers et al. 2011). By using this data, 
we roughly estimate the mass of dense molecular gas (Mdense) 
in the CND as follows. First, we convert L'^cn( 4 - 3 ) 
responding value of HCN(l-O), assuming HCN(4-3)/HCN(l- 
0) = 0.6±0.2 in the brightness temperature unit. This ratio is 
an averaged value over the CND of the nearby typical Seyfert- 
2 galaxy, NGC 1068, sampled at 100 pc resolution (Viti et al. 
2014). Note that NGC 1068 has a quite similar X-ray lumi¬ 
nosity and an Eddington ratio to NGC 7469. Then, we obtain 
^HCN(i-o) - (9.3±3.1) X 10^ K km s“^ pc^. Einally, assuming 
a conversion factor from L[^cn(i-o) ^dense to be 10 , which 
is derived for AGNs by Krips et al. (2008) based on IRAM 
30m observations, we estimate Mdense = (9.3±3.1) x 10^ Mq. 


The dynamical mass (Mdyn) of the CND, on the other hand, 
is calculated assuming a thin disk with Keplerian rotation as, 
Mdyn = 230 (r/pc) [v(r)^/(km s“^)] (sin i)~^ Mq . Here, v(r) is a 
rotation velocity at a radius r, and i is an inclination angle of 
the disk. Then, putting the EWHM ^ 200 km s“^ at the CND 
(1" = 330 pc diameter) and i = 45° (Davies et al. 2004), we 
obtain Mdyn = 7.6 x 10^ Mq. This is almost identical to the 
value derived in Eathi et al. (2015), who isolated the rotation 
velocity and the velocity dispersion based on their harmonic 
decomposition analysis of the observed velocity held. Our 
Mdyn is also consistent to those derived for r < 2.5" region, 
i.e., including the SB ring (6.5 x 10^ Mq; Davies et al. 2004). 

Taking them together, we can estimate the dense gas mass 
fraction in the CND of NGC 7469 as /gas = Mdense/Mdyn ^ 
12%. Recalling the fact that the above used mass conversion 
factor has a factor of^3 uncertainty (Krips et al. 2008), /gas 
eventually has the same level of uncertainty at this moment. 
Therefore, although the value of /gas ^ 12% is comparable to 
those estimated for low-luminosity AGNs (Sani et al. 2012) 
using the same conversion factor (i.e., 10 ), it is still difficult 
to speculate whether or not /gas has any relevance with the 
current mass accretion onto the central supermassive black 
hole. We need an accurate mass conversion factor estimated 
through multi-transitional analysis (i.e., estimating an accu¬ 
rate relative abundance of HCN with respect to H 2 ) to further 
constrain the /gas and its relevance to the nuclear activity. 

5. SPATIAL DISTRIBUTION OF THE DENSE GAS TRACERS 

We show in Eigures 3 to 5 the channel maps of the CO(3- 
2), HCN(4-3), and HCO'^(4-3) line emissions from the cen¬ 
tral 8" X 8" (2.6 kpc X 2.6 kpc) region of NGC 7469. We 
could not cover the whole CO(3-2) emission line as noted be¬ 
fore. Comparing with CO(2-l) map in Davies et al. (2004), 
we conclude that the emission from the south-west region of 
the map is missing. The HCN(4-3) and HCO‘^(4-3) emissions 
were detected (i.e., > 3 a) over velocity ranges of 4744-5165 
km s“^ and 4752-5150 km s“^ at the position-A, respectively, 
of which widths (EWZI ~ 400 km s“^) are almost consistent 
with each other. As for CO(3-2), the covered velocity range 
corresponds to ^ 70% of the total velocity range of HCN(4-3) 
and HCO'^(4-3) at the position-A. However, we can not see a 
clear turnover feature in the CO(3-2) spectrum at that position 
(Eigure 2). Therefore, we conclude that the EWZI of CO(3-2) 
emission line is signiflcantly wider than those of HCN(4-3) 
and HCO‘^(4-3) lines there. It might be possible that we de¬ 
tected CO(3-2) emission from an outer diffuse, high-velocity 
component around the CND. It might also be possible that the 
high velocity component of HCN(4-3) and HCO‘^(4-3) emis¬ 
sion lines are too faint to be signiflcantly detectable. On the 
other hand, the line widths at the knots of the SB ring are 
almost consistent among CO(3-2), HCN(4-3), HCO'^(4 - 3) 
emission lines (Eigure 2). The channel maps also show that 
the CND and the SB ring are bridged by a bar-like, or several 
spiral arm-like structures as morphologically discussed later. 
A detailed kinematic study of this region is presented in Eathi 
etal. (2015). 

We also show the maps of the CO(3-2), HCN(4-3), 
HCO‘^(4-3), and CS(7-6) velocity-integrated intensities in the 
central 2 kpc region of NGC 7469 in Eigure 6 . The peak po¬ 
sitions of the 860 jim continuum emission (the positions A-D 
in Eigure la) are again marked by the arrows in the CO(3- 
2) map. The dense gas tracers are also bright at these peak 
positions, indicating they are the sites of dust-obscured star 
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Table 4 

HCN(4-3) and HCO'''(4-3) line luminosities both in [Lq] unit and [K km s“* pc^] unit 


f-HCN(4-3) 

^HCO+(4-3) 

T' 

^HCN(4-3) 

T' 

^HCO+(4-3) 

[^q] 

[^0] 

[K km s“‘ pc^] 

[Kkms“‘ pc’] 

AGN position (7.95±0.34) x 10^ 

(7.01 ±0.27) X 10^ 

(5.58±0.24) X 10'’ 

(4.83±0.19) X 10'’ 

CND (2.10±0.09) X 10* 

(1.86±0.08) X 10* 

(1.47±0.06) X 10’ 

(1.28±0.06) X 10’ 


Note. — Row 1: Line luminosities within a synthesized beam (150 pc-resolution) placed at the position-A (i.e., AGN position) marked in Figure la. The 
luminosity was calculated following the formulae by Solomon & Vanden Bout (2005). Row 2: Line luminosities within the CND (i.e., the central ~ 1" = 330 pc 
diameter region). 



RA (J2000) 


Figure 3. Velocity channel maps of the CO(3-2) line emission in the central 8" x 8" (2.6 kpc x 2.6 kpc) region of NGC 7469. The central red cross in each 
channel indicates the position-A in Figure la (i.e., the AGN position). The velocity width of each channel is ~ 20 km s“^ and the central velocity of each channel 
is shown in the upper left corner. The synthesized beam (0.50" x 0.40", PA = 68.5°) is also plotted in the bottom left corner. Contour levels are -5, 3, 5, 10, 30, 
50, 70, 90, and 110 (t, where 1 a = 2.34 mJy beam“^. The negative contour is indicated by the dashed lines. Note that we could not observe channels with Vlsr 
< 4900 km s“^ due to our observational setup. Attenuation due to the primary beam pattern is corrected. 


formation. Two rather fainter knots are also seen in the CO(3- 
2) map, which are located northwest to the CND (marked as 
E and F). 

From these maps, it is apparent that the HCN(4-3), 
HCO‘^(4-3), and CS(7-6) emissions are significantly concen¬ 
trated towards the CND. A two-dimensional Gaussian fit to 
the HCN(4-3) image shows that the CND has an extent of 
(244.2±9.9) X (194.7±6.6)pc2 with PA = 58.8°±6.7°, which 
is similar to that of HCO'^(4-3) emission. This is the beam- 


convolved size, thus the actual CND must be smaller. We 
derived these uncertainties using Monte Carlo techniques by 
creating 10^ realizations of the image, each time adjusting 
the data by adding the normally-distributed (measured) noise 
randomly. Again using Figure 6, we estimate that the CND 
of NGC 7469 contains ^ 82 % and 53 % of the total emis¬ 
sion of HCN(4-3) and HCO'^(4-3) within the 18" field of view 
of ALMA (^6 kpc), respectively, suggesting their utility to 
probe the nuclear regions of external galaxies especially about 
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Figure 4. Velocity channel maps of the HCN(4-3) line emission in the central 8" x 8" (2.6 kpc x 2.6 kpc) region of NGC 7469. The central red cross in each 
channel indicates the position-A in Figure la (i.e., the AGN position). The velocity width of each channel is ~ 20 km s“^ and the central velocity of each channel 
is shown in the upper left corner. The synthesized beam (0.49" x 0.39", PA = 70.4°) is also plotted in the bottom left corner. Contour levels are -5, 3, 5, 10, 
20, and 30 a, where 1 cr = 0.55 mJy beam“^. Although we set the negative contours to be indicated by the dashed lines, they are not prominent in this hgure. 
Attenuation due to the primary beam pattern is corrected. 


HCN(4-3). 

Here, we expect that the missing flux is not signiflcant 
for the CND, although we recovered only 54±16 % of the 
HCO‘^(4-3) flux measured by APEX (Zhang et al. 2014) 
within the central 18" region, since its size 1") is much 
smaller than that of the maximum measurable angular scale 
of our observations, which was 7.1" (2.3 kpc). Since the 
apex’s HC0'^(4-3) spectrum of Zhang et al. (2014) is some¬ 
what noisy, we can expect the systematic flux uncertainty 
would be larger than the value mentioned in that paper (15% 
for absolute flux calibration uncertainty) additionally due to, 
e.g., baseline uncertainty. In other words, we expect the actual 
recovered rate can be higher than 54%. Unfortunately, we can 
not estimate the percentage of the recovered flux of HCN(4-3) 
and CS(7-6) since there are currently no signiflcant detection 
of these line emissions with single dish telescopes. 

In addition to the CND, we can see a partial circumnu- 
clear SB ring at a radius of 1.5"-2.5". Davies et al. (2004) 
also reported the existence of a bar or a pair of short, loosely 
wound spiral arms at a PA = 56° between the CND and the SB 


ring (connecting the positions B-A-D in Eigure 6a), based on 
their 0.7" CO(2-l) observations. Interestingly, that structure 
is only found in molecular emission; neither K-band (Gen- 
zel et al. 1995; Scoville et al. 2000) nor centimeter continuum 
(Condon et al. 1991 ; Orient! & Prieto 2010) image shows such 
a structure. This conflguration is quite unusual if the structure 
is a genuine molecular bar, since a corresponding stellar struc¬ 
ture should be easily visible in this scenario. Indeed, Davies 
et al. (2004) found no kinematic characteristics expected for 
a barred potential, i.e., an S-shaped contours in a isoveloc¬ 
ity held (e.g., Lundgren et al. 2004; Kohno et al. 1999) and 
a tilted-X shape in a position-velocity diagram along the (hy¬ 
pothesized) bar axis (e.g., Laine et al. 1999; Sakamoto et al. 
2000). Moreover, Eigures 4 and 5 reveal that there seem to 
be several connections between the CND and the SB ring in 
addition to the above hypothesized bar axis, e.g., channels at 
around 4850, 4980, and 5020 km s“^ Although we can not 
rule out the existence of a faint stellar bar that could not be 
detected in previous NIR observations, it might be possible 
that the nuclear region of NGC 7469 consists of several spiral 
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Figure 5. Velocity channel maps of the HCO'^(4-3) line emission in the central 8'' x 8'' (2.6 kpc x 2.6 kpc) region of NGC 7469. The central red cross in 
each channel indicates the position-A in Figure la (i.e., the AGN position). The velocity width of each channel is ~ 20 km s“^, and the central velocity of each 
channel is shown in the upper left corner. The synthesized beam (0.48" x 0.38", PA = 57.9°) is also plotted in the bottom left corner. Contour levels are -5, 3, 5, 
10, 20, and 30 cr, where 1 cr = 0.58 mJy beam“^. Although we set the negative contours to be indicated by the dashed lines, they are not prominent in this hgure. 
Attenuation due to the primary beam pattern is corrected. 


features like the case of, e.g., NGC 1097 (Prieto et al. 2005; 
Fathi et al. 2006; Davies et al. 2009). Since the lack of spatial 
resolution leads to blending several molecular clouds, which 
would result in an illusional bar, higher angular resolution and 
sensitivity observations of both CO and NIR emissions are 
necessary to settle this issue. 

6. SUMMARY AND CONCLUSION 

In this paper, we present high resolution (0.5" x 0.4"; 1" 
= 330 pc) ALMA band 7 (350 GHz band) observations of the 
submillimeter dense molecular gas tracers such as HCN(4-3), 
HCO'^(4-3), CS(7-6), and partially CO(3-2), in addition to the 
underlying 860 /im continuum emission, towards the central 
kpc region of the luminous type-1 Seyfert galaxy, NGC 7469. 
The region consists of the CND (central ^ 1") and the sur¬ 
rounding SB ring with a radius of L5"-2.5". We revealed the 
spatial distribution of these emissions. Thanks to the high spa¬ 
tial resolution, we can reliably measure the line fluxes at the 
CND and the SB ring separately. The obtained line ratios are 
directly compared with those in NGC 1097 (a low-luminosity 


type-1 Seyfert galaxy), to achieve insights of possible effects 
of AGN luminosity on the surrounding molecular material. 
The main results of our observations and conclusions are sum¬ 
marized as follows: 

- The spatial distribution of the 860 /im continuum emis¬ 
sion resembles well to those of centimeter and MIR 
continuum emissions, and dense molecular gas tracers. 
But their peak positions (including some knots in the 
SB ring) are shifted signiflcantly from the optical ones, 
suggesting the existence of dust obscured star forma¬ 
tion there. 

- The maximum position of the 860 /im continuum is 
identical to that of the VLA 8.4 GHz continuum, indi¬ 
cating it is the AGN position. The SED reveals that the 
860 /im continuum there is dominated by a thermal dust 
emission. We detected CO(3-2), HCN(4-3), HCO+(4- 
3), and CS(7-6) emissions at the same position. Other 
lines including the vibrationally excited HCN were not 
detected, although NGC 7469 is classifled as a LIRG. 
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Figure 6. Integrated intensity maps of (a) CO(3-2), (b) HCN(4-3), (c) HCO^(4-3), and (d) CS(7-6), in the central ~ 2 kpc region of NGC 7469. The rms noises 
(1 fj) are 0.20, 0.06, 0.06, and 0.06 Jy beam“^ km s“^ in (a), (b), (c), and (d), respectively. Note that we could not observe CO(3-2) emission at Vlsr < 4900 km 
s“^ due to our observational setting, thus the displayed map is incomplete (south-east part is deficient). Some bright molecular knots are marked as A-F in the 
CO(3-2) map (the positions A-D are the same as the ones marked in Figure la). The white filled ellipses indicate the synthesized beams (~ 0.50" x 0.40"). The 
central crosses indicate the AGN position. Contours are: (a) -5, 3, 5, 10, 20, 30, ..., and 130 a, (b) -5, 3, 5, 10, 20, 30, ..., and 70 a, (c) -5, 3, 5, 10, 20, 30, 40, 
and 50 a, (d) -5, 3, and 5 a, respectively. Negative contours are indicated by dashed lines. See also Section 2 for the velocity ranges integrated over. 


- The high Rucn /hco+ Rucn /cs are only found in the 

CND, which is consistent to the findings in 113, i.e., 
AGNs tend to show higher ^hcn/hco+ 1) and/or 
Rucn/cs 10) than those in SB galaxies (submm- 
HCN enhancement). 

- Although the actual physical mechanisms to enhance 
7?hcn/hco+ in AGNs are still unclear, we found this ratio 
is significantly lower in NGC 7469 (1.11 ±0.06) than in 
NGC 1097 (2.0±0.2) despite the more than two orders 
of magnitude higher X-ray luminosity of NGC 7469 
than that of NGC 1097. The close vicinity of the lumi¬ 
nous AGN of NGC 1068 also shows a comparable ratio 

1.5) to NGC 7469. Some other heating mechanisms 
than X-ray (e.g., mechanical heating) would contribute 
significantly for shaping the chemical composition in 
NGC 1097. 

- The CO(3-2) emission, which is not fully observed 
though, is distributed over both the CND and the SB 
ring, whereas the HCN(4-3) and the HCO'^(4-3) emis¬ 
sions are significantly concentrated toward the CND 

82% and 53% of the total fiux inside the 18" field 
of view of ALMA, for HCN and HCO”^, respectively), 
suggesting their utility to probe the nuclear regions of 
galaxies. 


- Between the CND and the SB ring, we found several 
spiral-like structures, which would be transporting gas 
from the SB ring to the CND. This might be a more 
preferable scenario than a bar. We should test this ten¬ 
tative view by further high resolution observations and 
kinematic analysis. 

From this work, we achieved supportive evidence for the 
submm-HCN enhancement in AGNs tentatively proposed by 
113. We will further investigate this trend by increasing the 
number of sample galaxies, and extensively discuss the pos¬ 
sible cause of the HCN-enhancement in our succeeding paper 
(Izumi et al. submitted). 

We thank the anonymous referee for very kind and 
helpful comments for improving the paper. This 
paper makes use of the following ALMA data: 
ADS/JAO.ALMA#2012.L00165.S. ALMA is a partner¬ 
ship of ESO (representing its member states), NSF (USA), 
and NINS (Japan), together with NRC (Canada) and NSC 
and ASIA A (Taiwan), in cooperation with the Republic of 
Chile. The Joint ALMA Observatory is operated by ESO, 
AUI/NRAO, and NAOJ. The National Radio Astronomy 
Observatory is a facility of the National Science Founda¬ 
tion operated under cooperative agreement by Associated 
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Universities, Inc. We used data based on observations with 
the NASA/ESA Hubble Space Telescope, and obtained from 
the Hubble Legacy Archive, which is a collaboration be¬ 
tween the Space Telescope Science Institute (STScI/NASA), 
the Space Telescope European Coordinating Eacility (ST- 
ECE/ESA) and the Canadian Astronomy Data Centre 
(CADC/NRC/CSA). In addition, this research has made use 
of the NASA/IPAC Extragalactic Database (NED) which 
is operated by the Jet Propulsion Laboratory, California 
Institute of Technology, under contract with the National 
Aeronautics and Space Administration. T. 1. was supported 
by the ALMA Japan Research Crant of NAOJ Chile Observa¬ 
tory, NAOJ-ALMA-0029. T. 1. and H. U. are thankful for the 
fellowship received from the Japan Society for the Promotion 
of Science (JSPS). 
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Figure 7. The heating rate per hydrogen atom due to X-ray absorption (red solid) and to soft-UV absorption (blue dashed) as a function of radius from the 
center, under the ID slab geometry. We assume that a Mbh = 10^ Mq black hole is emitting 3% of its Eddington luminosity at 1-10 keV. The soft-UV radiation 
field is taken as Go = 10^-^ (Habing unit). The adopted gas density is nu = 10^ cm“^. In this configuration, the estimated size of the XDR is r ~ 42 pc. 

APPENDIX 

A. THE SIZE OF AN XDR 

Here we estimate the size of an XDR, which is defined as the largest distance from an AGN where X-ray heating due to the 
AGN dominates over UV heating due to SB, by employing a toy model. 

For the X-ray incident flux, we adopt a spectral shape of the form, 

/^(£) = Fo(y^) exp(-£/£,), (Al) 

where Fb is a scaling constant, E = hu eV, ct is a power-law index, Ec is a cut-off energy. We set ct = 1.0 (e.g., Madejski et al. 
1995; Peterson 1997) and Ec = 100 keV (e.g., Madejski et al. 1995). We hereafter consider the 1-10 keV photons for simplicity. 
The softer X-ray can be easily absorbed at the edge of the molecular cloud and thus neglected here. Considering the situation 
that an AGN with r = 0.1 pc and Mbh = 10^ M©, is emitting isotropically at 3% of its Eddington luminosity at this 1-10 keV 
range, we scaled Eq. Subsequent radiative transfer is calculated by adopting a ID slab geometry with a uniform gas density of 
= 10^ cm“^. Note that the AGN of NGC 7469 hosts a ^ 10^ Mq black hole with the Eddington ratio of ~ 0.3 (see also Section 
1). Thus, our estimation corresponds to the case where 10% of the bolometric AGN luminosity of NGC 7469 is going into the 
1-10 keV band. Then, an energy deposition rate per particle (E/x; Maloney et al. 1996) is 

/•lOkeV 

/7x= / apa(^)F(^)exp(-r)J£, (A2) 

7lkeV 

with an optical depth of 

r = apa(^).AH(r). (A3) 

Here, we use the photoelectric absorption cross section crpa(£’) of Morrison & McCammon (1983), which is estimated for 1-10 
keV range. The attenuating hydrogen column density from the AGN to the point at the distance of r pc is indicated by Ar. 
Assuming the heating efficiency of X-ray photons to be 30% (e.g., Maloney et al. 1996; Meijerink & Spaans 2005), an X-ray 
heating rate is calculated as a function of a radius from the AGN (Figure 7). On the other hand, we adopt a typical energy of 10 
eV and a cross section of 2.78 x 10“^^ cm^ for soft-UV photons from the SB activity (Meijerink & Spaans 2005). This SB is 
assumed to produce a soft-UV radiation field of Go = 10^'^ (Habing unit). Note that Rosenberg et al. (2014a,b) used this Go to 
reproduce /upper ^ 4 emission lines (i.e., lines tracing cold molecular gas) observed in NGC 253 and Arp 299. Assuming 
that the UV heating efficiency is ^ 0.3%, the heating rate due to SB-induced UV radiation is also shown in Figure 7. 

Although the above estimation strongly depends on the adopted parameters and the value of Go is not well constrained in the 
circumnuclear region of Seyfert galaxies at this moment, we can roughly estimate that the size of the XDR from Figure 7, which 
is r ~ 42 pc in this case. This size is well smaller than the spatial resolution of our ALMA observations (150 pc). 

B. ATTENUATED X-RAY FLUX 

Using the Equations (Al) and (A3), we can calculate an attenuated 1-10 keV X-ray flux as a function of radius. The cases of 
NGC 7469 and NGC 1097 are investigated. Again we used the photoelectric absorption cross section calculated by Morrison & 
McCammon (1983). We adopt the same values as indicated in Appendix- A for this calculation, except for Eq. We this time scaled 
Eq by using the 2-10 keV luminosity of NGC 7469 and NGC 1097 (e.g., Liu et al. 2014), a = 1.0, and assuming their AGN sizes 
to be r = 0.1 pc. Then, the attenuated X-ray fiux is shown in Figure 8. 
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Figure 8. ID calculation of the 1-10 keV X-ray flux as a function of radius from the AGN. Attenuation due to obscuring material (Morrison & McCammon 
1983) is taken into account. The cases of NGC 7469 (red solid) and NGC 1097 (blue dashed) are shown. The horizontal dot-dashed line indicates the flux of 10 
erg s“^ cm“^, which is a required value by Meijerink et al. (2007) to realize f?HCN/HCO+ > 1- Thus, the largest radius that can keep Ti_iokeV ^ 10 is ~ 35 pc and 
~ 18 pc, for NGC 7469 and NGC 1097, respectively. 






